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Cooling of a Bose-Einstein Condensate by spin distillation 
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We propose and experimentally demonstrate a new cooling mechanism leading to purification of 
a spinor Bose-Einstein Condensate (BEC). Our scheme starts with a BEC polarized in the lowest 
energy spin state. Spin excited states are thermally populated by lowering the single particle energy 
gap set by the magnetic field. Then these spin-excited thermal components are filtered out, which 
leads to an increase of the BEC fraction. We experimentally demonstrate such cooling for a spin 
3 ®^Cr dipolar BEC. Our scheme should be applicable to Na or Rb, with perspective to reach 
temperatures below 1 nK. 

PACS numbers: 37.10.De, 03.75.Mn,75.30.Sg 


In the last three decades, laser cooling and evapora¬ 
tive cooling have led to major advances in atomic and 
molecular physics, in particular in the fields of precision 
measurements, atomic clocks, and quantum degenerate 
gases [lilll- Nowadays, the hope to stucW magnetic cor¬ 
relations of atoms in optical lattices d, Q , and the pos¬ 
sible connections to exotic superconductivity are major 
motivations to obtain systems with lower entropies than 
currently available Q. It is therefore important to find 
new ways to remove entropy in degenerate quantum Bose 
or Fermi gases loaded in optical lattices [6l-l9|. 


In typical lattice experiments, atoms are first cooled 
by evaporative cooling, and then loaded in the periodic 
potential. Unfortunately, evaporative cooling ceases to 
be efficient when the temperature is significantly smaller 
than the interaction ener gy, which currently sets an ul¬ 
timate limit for entropy [lOj. A number of alternative 
cooling schemes have been suggested and studied 11 M 
but up to now the best phase space densities are still due 
to evaporative cooling. Here, we propose to use the spin 
degrees of freedom to efficiently store and remove entropy 
in partially Bose-condensed gases to reach temperatures 
below the current limitations set by evaporative cooling. 


Our proposal starts with a sample polarized in the low¬ 
est energy spin state. We engineer a thermodynamic cy¬ 
cle (see fig. H^)) wherein the magnetic field is first low¬ 
ered to trigger depolarization of the gas, and the resultant 
spin-excited states are then filtered out of the trap. As 
shown in fig. [It), this cycle reproduces a BEC polar¬ 
ized in the lowest energy state with an increased conden¬ 
sate fraction, provided the initial thermal fraction is low 
enough. The gain in phase space density directly follows 
from Bose statistics: as the condensate forms in the low¬ 
est energy single-particle spin state [l7|, spin filtering of 
the excited spin states obtained after depolarization in¬ 
troduces a loss which is specific to thermal atoms. At low 
magnetic field, spin filtering typically leads to a decrease 
of the thermal fraction, and hence of entropy, by a factor 



FIG. 1: a) Principle of the proposal, b) Density profiles after depo¬ 
larization and Stern Gerlach procedure showing a ®^Cr BEC forming 
only in rris — —3. Circles are experimental data, the solid line is result 
of bimodal fit, and the dashed line only fits thermal fractions, c) Pu¬ 
rification of a ®^Cr BEC. After depolarization occurring at a B field of 
1 mC, we measure the final condensate fraction and hence the effective 
final reduced temperature (red squares), as a function of the initial 
reduced temperature. Error bars show statistical uncertainties. The 
(red) solid line is the result of our model. The (black) straight solid 
line (of unity slope, crossing the origin) is a reference (i.e. no cooling). 
The inset illustrates cooling by successive cycles. 


~ 2. The cycle can then be repeated. We find that this 
cooling efficiency has no fundamental limit other than 
that associated with technical noise. 

To demonstrate this new cooling strategy, we consider 
two specific systems. The first system is a chromium 
dipolar gas, which is sensitive to the linear Zeeman effect, 
because magnetization is free [l^ . Bose-Einstein conden¬ 
sation occurs in the lowest energy spin state rris = —3 
(above a critical B field, see below). Then, thermal pop¬ 
ulation of TOs > —3 spin-excited states due to dipole- 
dipole interactions (see fig. [T]d)) can be used for cooling. 
The cooling limit arises from technical difficulties at con¬ 
trolling the magnetic field at the 100 ^G level. Despite 
these difficulties, this paper provides a proof-of-principle 
experiment, demonstrating the efficiency of spin filtering 
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to purify a chromium BEC in an optical dipole trap. 

The second system is a spinor s = 1 BEC such as Rb 
or Na, prepared in rUg = 0. Spin-changing collisions as¬ 
sociated with the difference in scattering lengths in the 
molecular potentials S = 0 and S = 2 redistribute popu¬ 
lation between rus = 0 and iris = ±1 at constant magne¬ 
tization. The (positive) quadratic Zeeman effect provides 
an energy shift q between a pair of atoms in the ms = 0 
state and a pair of atoms in states mg = ± 1 , which fa¬ 
vors BEC in mg = 0. Thermal population of spin-excited 
states is possible as long as fc^T > q, which sets a prac¬ 
tical limit to cooling in the pK range. This limit ensures 
that spin fluctuations of mesoscopic polar BECs are neg¬ 
ligible 0 . 

Before discussing our experimental results, we now 
turn to our theoretical model. The suggested experi¬ 
ments start with a finite temperature polarized BEC in 
the lowest energy spin state to™™. The magnetic field 
is then rapidly (diabatically) reduced to allow depolar¬ 
ization. To compute the equilibrium state after depolar¬ 
ization, we assume thermodynamic equilibrium (insured 
by collisions), and conservation of the total atom num¬ 
ber Ntot and total energy. We use the thermodynam¬ 
ics of non-interacting atoms, assuming that mean-field 
interactions can be neglected when evaluating the equi¬ 
librium state, so that spin-dependent Bose occupation 
factors read: 


condensate 

fractions 


a) 


S,/S, 


b) 


fk,m^ (m) — 


exp [P {ek,Tn, - m)] - 1 


( 1 ) 


where p, is the chemical potential, /3 = l/ksT, and ek,ms 
is the single-particle energy of the (trapped) states la¬ 
beled by the index fc, in the Zeeman state TOg. When 
magnetization is free (case of dipolar particles) ek,ms in¬ 
cludes the linear Zeeman effect, whereas when magneti¬ 
zation is fixed (as for example for Rb and Na atoms), the 
linear Zeeman effect is gauged out and ek,ms only includes 
the quadratic Zeeman effect q. In both cases, we spec¬ 
ify that the lowest energy spin state to™™ (to™™ = —3 
for Cr with free magnetization, to™™ = 0 for Rb or Na 
at zero magnetization) has = 0 for the lowest 

trapped states k ^ 0 (which sets the otherwise arbitrary 
energy scale). We assume that a BEC is present in state 
TO™™, so that fj, = 0. 

In the case of a 3D harmonic trap, and in the thermo¬ 
dynamic limit, the thermal population in each spin state 
reads: 


fVm, ='^fk,mS0) « 93 (e 


( 2 ) 


with gn the polylogarithm function of order n, and ui the 
geometric average of the angular trapping frequencies of 
the 3D harmonic trap. Additional atom are condensed 
in state to™™. 



Condensate fractions for fi — 0.7 after depolarization, with (solid, or¬ 
ange) or without (dashed, red) spilling, as a function of By. b) Ratio of 
final (after spilling) and initial mean entropy per atom as a function of 
initial temperature, for By=0.3 mG (dashed) and By=0.1 mG (solid). 
The horizontal dotted line corresponds to the maximal gain, reached 
at By=0 and Ti — 0 (given by eq. The vertical lines indicate the 

qualitative limit for cooling, set by — ksT. The dots (black) 

correspond to experimental data taken at By = 0.2±0.1 mG, including 
statistical uncertainties. 

The total energy of the system J2k,m, fk,ms (0)efc,m„ 
is the sum of trap (kinetic and potential) and magnetic 
energies: 


Etrap ~ ^3fcBTg4(e 
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The final temperature Tf and BEC atom number Ncj 
are thus derived, and can be compared to initial values Ti 
and Nc^i- Any atom in spin-excited state TOg ^ to™™ can 
be removed (by means of magnetic field gradients, micro- 
wave transitions, and/or a resonant push laser beam). 
We define two post-depolarization BEC fractions, with or 
without spin filtering the excited states, which are respec¬ 
tively /2 = Ncj/{Ncj + N^^ir^j), and /i = Ncj/Ntot 

{Ntot = -|- N^min i). 

As shown in fig. [It), we find that this procedure should 
lead to an increase in condensate fraction /2 only when 
the initial condensate fraction fi = l — {Ti/Tdf' is large 
enough. We interpret this complex behaviour as the con¬ 
sequence of the competition between two effects, (i) As 
population in spin-excited states is purely thermal, spin 
filtering leads to purification of the BEC. (ii) As the gas 
depolarizes and the number of thermal atoms in to™™ 
decreases, the BEC must melt to maintain saturation. 
This competition can lead to an increase in condensate 
fraction, because BEC atoms have zero energy. Hence, 
melting of the BEC cools the thermal gas in to™™, which 
can then be saturated at a lower temperature, as already 
observed in and [ 2 lj |. 

At low enough T, the cooling process becomes most 
efficient for B^- 0 (9 0 ), and in that limit we obtain 

the following simple expression: 


1-/2 = 


1 -/. 


(2s + 1)3/4 _2s(l_y) 


( 4 ) 
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Magnetic field (mG) 


Magnetic field (mG) 


Fldj. Si (Color online) Experimental results for spin spilling with 
the s — 3 ^^Cr BEC, for an initial BEC fraction of 0.65, and different 
values of B/. a) Number of atoms after spin spilling; in the BEC (top), 
in the thermal component, b) BEC fraction deduced from a). The solid 
lines are guides for the eye. 


Therefore, at zero B field, while /i = 1 — (2s + — 

fi) < h > ft for large enough /, (if fi > 0.45 for 
s = 3). 

It is interesting to assess the efficiency of this cooling 
procedure. The best way to do so is to consider the re¬ 
duction of entropy after each cooling cycle. Below Tc, 
all entropy is contained in the thermal fraction, and each 
thermal atom has a temperature-independent entropy of 
3.6 fcs Q. As a consequence, selectively removing ther¬ 
mal atoms is an excellent way to remove entropy. For 
B=0, ^ ^ (2s-f-i)3/4 when fi ^ 1, with Sij the 

mean entropy per atom before depolarization, and after 
spilling, respectively. In figElD) the entropic efficiency of 
one cooling cycle is shown for different values of the B 
field, as a function of Ti, in the case of s = 3 Cr. ^ 
remains significantly smaller than 1 even at very low Ti, 
provided a low enough B field: a reduction of entropy by 
a factor of at least two is obtained until ksTi ~ 

This corresponds to Ti ~ 0.1 x Td = 30 nK in our exper¬ 
iment. With s = 1 Na or Rb, this limit can be pushed 
to even lower Ti. For example for Na, with an energy 
difference of 280 Hz/G^ between {rus = 0;ms = 0} and 
{nis = —l;ms = 1}, temperatures in the 100 pK range 
can be achieved with B fields below ~ 100 mG. 

We have also estimated the effect of interactions, by 
including in ek,ms the effect of spin-dependent contact 
interactions within the Bogoliubov approximation (^ . 
We find that interactions between particles do not modify 
the general picture described above, as long as the BEC 
remains polarized in a well defined Zeeman state. For 
the case of Na, spin-dependent interactions favor a polar 
TOs = 0 BEC. For Rb, at typical densities of 10^"^ at.cm“^, 
a quadratic energy shift of typically 10 Hz, provided by a 
magnetic field of 380 mG, is sufficient to insure that the 
BEC is polar. For the case of chromium (Rb), the BEC is 
polarized in rris = —3 (rUs = 0) above a critical magnetic 


field Be set by spin-dependent interactions [2J, [2^. In 
the experiment reported here. Be ~ 100 /rC. 

Finally, we also find a subtle effect associated with 
demagnetization cooling (a process previously demon¬ 
strated in the group of T. Pfau in the thermal regime 
[l5|l. When the experiment is performed at non-zero 
magnetic field, the cooling associated with the exchange 
between magnetic and kinetic energy can lead to an in¬ 
crease in the BEC atom number even without spin filter¬ 
ing: Nej > Nc,i, and hence fi > fi, as shown in fig[2ti). 
However, this effect is small compared to the two major 
effects described above (and smaller than our experimen¬ 
tal error bars). The increase in BEC atom number occurs 
despite the necessary increase in total entropy, because 
spin-excited states have more entropy per particle than 
the saturated (thermal) lowest energy state. This bears 
similarity with the reversible production of a BEC when 
tuning the anharmonicity of a trap 16l |. 

We now present our experimental results demonstrat¬ 
ing purification of a ®^Cr BEC. Condensation is obtained 
in the absolute ground state = —3, by evaporative 
cooling in a crossed optical dipole trap obtained from a 
100 W, 1075 nm, IR laser 26[. The final trap frequencies 
at the end of evaporation are varied in order to change 
the initial BEC fraction fi at the beginning of the spilling 
process. For example, for ujx,y,z = 2 tt x (250,300, 215) 
Hz (measured through parametric excitation, with 5% 
uncertainty), we find fi = 0.5. This is in agreement with 
a critical temperature for BEC Tc = 310 nK for 2.10'’’ 
atoms, and a measured temperature of 250 nK. 

The magnetic field, which is initially Bi=40 mC, is 
then lowered in 50 ms to a final value Bj, low enough to 
trigger depolarization {ggs^f — ksT). To calibrate Bf, 
we first search for the B field where maximal depolariza¬ 
tion occurs, by varying all three B components with steps 
of 100 fj,G, which defines the zero B field. We then leave 
two components of the B field unchanged, and vary one 
component, which is calibrated by use of radio-frequency 
spectroscopy. We then let the cloud evolve at B/ for 150 
ms. This is long enough for inelastic dipolar collisions to 
ensure depolarization, as the typical rate for these colli¬ 
sions is 15 s“^ for a magnetic field in the 0-1 mC range 


27|. 


The spin filtering procedure consists of i) raising the 
magnetic field to 30 mC; ii) changing the trap in or¬ 
der to spill out atoms with rUs > —3. For that we 
decrease the IR power and add a vertical B field gra¬ 
dient. In this trap, mg = ~3 atoms remain trapped, 
while TOs > — 3 atoms fall, because optical and magnetic 
forces are insufficient to compensate gravity for these 
atoms. We checked through Stern-Gerlach analysis that 
indeed the losses induced by the new trap are negligible 
for rUs = —3, but almost total for rUg > —3. In addition, 
we checked that no evaporation is induced by the filter¬ 
ing procedure. The spin filtering is performed rapidly 
enough (in 60 ms) to ensure that no dipolar relaxation 
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fraction after spilling, / 2 , is plotted as a function of B/. Lines are 
predictions of our non-interacting model. For large B/, /2 reaches fi. 

occurs. Finally, we recover the initial (purely) optical 
trap at and measure the new BEC fraction, / 2 . 

To accurately measure BEC fractions, we release the 
atoms from the trap and take an absorption picture af¬ 
ter a time of flight of 5 ms. We perform a two-stage 
analysis of the images. A first fit is used to measure the 
total number of atoms, and deduce the critical tempera¬ 
ture Tc- A second bi-modal fit is then performed, where 
the temperature is a free parameter which sets both the 
width of the gaussian describing thermal atoms, and the 
BEC fraction following non-interacting Bose thermody¬ 
namics predictions. This procedures ensures measure¬ 
ment of thermal fractions down to a few percent. Results 
are shown in fig l3] for an initial temperature (correspond¬ 
ing to fi = 0.65) low enough that the spilling is efficient. 
At the lowest B/, the spilling procedure induces a rela¬ 
tive loss in thermal atoms larger than that for condensed 
atoms. This figure illustrates the competition between 
the two effects described above, which results here in an 
increase of the BEC fraction / 2 . 

Figure |4] shows our complete set of data, for different 
initial condensate fractions fi. We observe different be¬ 
haviors for /2 at the lowest B/. For the smallest fi, /2 
gets smaller than fi. On the other hand, when fi is large 
enough, /2 gets significantly larger than fi. All these 
results agree with our simple non-interacting model. In 
terms of entropy, the measured efficiency of the scheme 
(see fig[2|) also agrees with our model; in particular, the 
reduction of entropy is more pronounced as Ti decreases. 

In conclusion, we propose an efficient cooling mecha¬ 
nism using on the spin degree of freedom. It is based 


on redistribution of entropy among the spin states, while 
recent theoretical proposals or experiments rely on 
spatial redistribution [7|, as is for example the case for 
spin gradient demagnetization cooling [8|, or distillation 
of singly-occupied sites in optical lattices [9| . In addition, 
cooling by spin filtering can be repeated an arbitrary 
number of times. Since each cycle leads to typically a 
factor of two reduction in mean entropy per atom, we 
foresee that this scheme could indeed be a way to reach 
new regimes of deep degeneracy. As very pure conden¬ 
sates are obtained, measuring small thermal fractions will 
be a challenge. It will then be advantageous to use the 
spin degrees of freedom for thermometry, as explored in 
[isj l , and discussed in : counting atoms in spin-excited 
states offers a background-free measurement, contrary to 
momentum distributions measurements. Finally, one of 
the interesting pending question is whether this scheme 
will help to remove entropy for BECs loaded in optical 
lattices, e.g. in the superfluid shells surrounding the Mott 
plateaux characteristic of the typical wedding cake dis¬ 
tribution [ 2 ^ . 
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